Kondo decoherence: finding the right spin model for iron impurities in gold and silver 
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We exploit the decoherence of electrons due to magnetic impurities, studied via weak localization, to resolve 
a longstanding question concerning the classic Kondo systems of Fe impurities in the noble metals gold and 
silver: which Kondo-type model yields a realistic description of the relevant multiple bands, spin and orbital 
degrees of freedom? Previous studies suggest a fully screened spin S Kondo model, but the value of S remained 
ambiguous. We perform density functional theory calculations that suggest S — 3/2. We also compare previous 
and new measurements of both the resistivity and decoherence rate in quasi 1 -dimensional wires to numerical 
renormalization group predictions for S = 1/2, 1 and 3/2, finding excellent agreement for 5 = 3/2. 

PACS numbers: 73.23.-b, 75.20.Hr, 72.70.+m, 73.20.Fz 



Introduction. — The Kondo effect of magnetic impurities in 
non-magnetic metals, e.g. Mn, Fe or Co in Cu, Ag or Au, 
first manifested itself in the early 1930's as an anomalous rise 
in resistivity with decreasing temperature, leading to a resis- 
tivity minimum 0]]. In 1964 Kondo explained this effect J2l 
as resulting from an antiferromagnetic exchange coupling be- 
tween the spins of localized magnetic impurities and delocal- 
ized conduction electrons. 

However, for many dilute magnetic alloys a fundamental 
question has remained unresolved to this day: which effec- 
tive low-energy Kondo-type model yields a realistic descrip- 
tion of the relevant multiple bands, spin and orbital degrees 
of freedom [3[]? Cases in point are Fe impurities in Au and 
Ag, the former being the very first magnetic alloy known to 
exhibit an anomalous resistivity minimum |lj]. Previous at- 
tempts to fit experimental data on, for example, Fe impurities 
in Ag (abbreviated as AgFe) with exact theoretical results for 
thermodynamics, by assuming a fully screened low-energy ef- 
fective Kondo model flBt], have been inconclusive: specific 
heat data is absent and the local susceptibility of Fe in Ag 
obtained from Mossbauer spectroscopy Q6Q indicated a spin of 
S = 3/2 while a fully screened S = 2-model has been used to 
fit the temperature dependence of the local susceptibility 10]. 

A promising alternative route to identify the model for Fe 
in Au or Ag is offered by studying transport properties of high 
purity quasi-one dimensional mesoscopic wires of Au and Ag, 
doped with a carefully controlled number of Fe impurities by 
means of ion-implantation 



decoherence. By checking model predictions for both ef- 
fects against experimental observations over several decades 
in temperature, decoherence can thus be harnessed as a highly 
sensitive probe of the actual form of the effective exchange 
coupling. Experiments along these lines III U Il2ll were con- 



sistent with a Kondo model in which the impurity spin is fully 
screened and inconsistent with underscreened or overscreened 



Kondo models III 111 . A consistent description of both resistiv- 
ity and decoherence measurements using the simplest fully 
screened Kondo model, the S — 1/2 single-channel Kondo 
model, was, however, not possible: different Kondo scales 
were required for fitting the resistivity and decoherence rates 
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12L 113D. Magnetic 
impurities affect these in two different ways. Besides caus- 
ing the afore-mentioned resistivity anomaly, they also make 
an anomalous contribution 7 m (T 1 ) to the electronic phase de- 
coherence rate 7^ (T) measured in weak (anti)localization: an 
itinerant electron which spin-flip-scatters off a magnetic im- 
purity leaves a mark in the environment and thereby suffers 



In this Letter we address the above problem via the follow- 
ing strategy: (i) We carry out density functional theory cal- 
culations within the local density approximation (LDA) for 
Fe in Au and Ag to obtain information that allows us to pre- 
scribe a low-energy effective model featuring 3 bands cou- 
pling to impurities with spin S = 3/2. (ii) We calculate the 
resistivity /o m (T) and decoherence rate j m (T) due to mag- 
netic impurities for three fully screened Kondo models, with 
n = 2S = 1, 2 and 3, using Wilson's numerical renor- 
malization group (NRG) approach, (iii) We compare these 
predictions to experimental data: extracting the characteristic 
Kondo temperature Tjf for each choice of n from fits to p m (T) 
and using these Tjf to obtain parameter-free predictions for 
7 m (T), we find that the latter agree best with experiment for 
n = 3. 

LDA calculations. — Fully relaxed densit y fu nctional the- 
ory calculations employing the VASP code 11411 showed that 
low-symmetry Fe configurations (split-interstitials 1 130 are 
energetically unfavorable: Fe impurities prefer an environ- 
ment with cubic symmetry. As the calculated defect formation 
energy of an Fe interstitial was found to be about 2 eV higher 
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FIG. 1: (color online). The d-level local density of states (LDOS) of 
substitutional Fe in Ag and Au within spin-polarized LDA (a,b) and 
LDA+U (c,d), with inclusion of spin-orbit interactions, and show- 
ing the e s (red) and tag (black) components of the d-level LDOS of 
FeAg (left panels) and FeAu (right panels). Majority/minority con- 
tributions are shown positive/negative. Legends give the spin ([is) 
and orbital (pl) magnetic moments in units of the Bohr magneton 
Pb, and the splitting (A) between the e g and t2 g components of the 
d-level LDOS. 



than the energy of a substitutional defect, we discuss the lat- 
ter case in the following. This is in line with experiments 
on Fe-implantation in AgAu alloys, where only substitutional 
Fe-atoms are found 111 611 . 

Fig. Q] shows the d-level local density of states (LDOS) of 
substitutional Fe in Ag and Au, obtained by spin-polarized 
calculations using a 108 atom supercell, with similar results 
being found for a 256 atom supercell. The cubic local sym- 
metry leads to e g (doublet) and t2 g (triplet) components with 
a e g -t 2g splitting, A > 0.15 eV in LDA (Fig. Ea-b)). The 
widths T c and T t2 of the e g and t2 g states close to the Fermi 
level (Ep) are of the order of 1 eV, resulting from a substantial 
coupling to the conduction electrons. The large t2 g component 
at E-p persists within LDA+U (Fig.[TJc-d) using U = 3 eV and 
a Hund's coupling Jh = 0.8 eV). 

The spin and orbital moments are given in the legends of 
Fig. Q] (spin-polarized Korringa-Kohn-Rostoker calculations 
yielded similar values lfl6Tl ): Within spin-polarized LDA a 
large spin moment fig of approximately 3-3.1 fiB forms spon- 
taneously, consistent with Mosbauer measurements that give 
3. 1-3.2 fiB for the spin moment for Fe in Ag Q|. In contrast, 
there is no tendency for a sizable orbital moment (or a Jahn- 
Teller distortion). The small orbital moments fiL of < 0.1 fiB 
(consistent with experimental results IU7I1 ) arise only due to 
weak spin-orbit coupling. We therefore conclude that the or- 
bital degree of freedom is quenched on an energy scale set by 
the width T t2s of the ta g orbitals. Moreover, since the spin- 
orbit splitting of the localized spin in the cubic environment 
is proportional to p\, it is tiny, well below our numerical pre- 
cision of 0.01 me V, and, therefore, smaller than the relevant 
Kondo temperatures. 



Low-energy effective models. — The above results justify 
formulating an effective low-energy model in terms of the 
spin-degree of freedom only. The large spin moment fig of 
3-3.1 fiB suggests an effective spin S = 3/2. Our LDA results 
thus imply as effective model a spin-3/2 3-channel Kondo 
model, involving local and band electrons of t2 g symmetry. 
An alternative possibility, partially supported by the large (al- 
most itinerant) t 2g component at Ep, would be to model the 
system as a spin-1 localized in the e g orbitals, that is perfectly 
screened by two conduction electron channels of e g symmetry. 
This spin is then also coupled to (almost itinerant) t2 g degrees 
of freedom via the ferromagnetic Jh- At high temperature, 
the latter binds an itinerant i 2g spin 1/2 to the local spin 1 to 
yield an effective spin-3/2, consistent with the spin-moment 
of 3-3.1 fiB obtained within LDA, whereas in the low temper- 
ature limit, the irrelevance of Jh under renormalization 
leads to the stated effective spin-1, 2-band model. Though 
such a model is well justified only for Jh <C r t2g , which 
is not the case here where Jh ~ T t2 , our LDA results do 
not completely exclude such a model. To identify which of 
the models is most appropriate, we shall confront their pre- 
dictions with experimental data below. 

We thus describe Fe in Ag and Au using the following fully 
screened Kondo model: 



yCkc 



(1) 



kaa 



It describes n channels of conduction electrons with wave vec- 
tor k, spin a and channel index a, whose spin density Ti a s a at 
the impurity site is coupled antiferromagnetically to an Fe im- 
purity with spin S = n/2. Whereas our LDA results suggest 
n = 3, we shall also consider the cases n = 1 and 2. 

NRG calculations. — The resistivity p m (T) and decoher- 
ence rate 7 m (T) induced by magnetic impurities can be ob- 
tained from the temperature and frequency dependence of the 
impurity spectral density 1118, Il9ll. We have calculated these 
quantities using the NRG 12 01 l2ll 12211 . While such calcula- 
tions are routine for n = 1 and 2 12111 . they are challenging 
for n = 3. Exploiting recent advances in the NRG ll20Tl we 
were able to obtain accurate results also for n = 3 (using a 
discretization parameter of A = 2 and retaining 4500 states 
per NRG iteration). 

Fig. |2] shows p m (T) and 7 m (T) for n = 2S = 1,2 and 
3. For T > Tj|, enhanced spin-flip scattering causes both 
p m (T) and j m (T) to increase with decreasing temperature. 
For T < T-H the effective exchange coupling becomes so 
strong that the impurity spins are fully screened by conduc- 
tion electrons, forming spin singlets, causing p m (T) to sat- 
urate to a constant and^m(T) to dro p to zero. While these 



effects are well-known IbSS IS 11, 12], it is of central im 



portance for the present study that they depend quite signifi- 
cantly on 5 = n/2, in such a way that conduction electrons 
are scattered and decohered more strongly the larger the local 
spin S: With increasing S, (i) both resistivities and decoher- 
ence rates decay more slowly with T at large temperatures 
Tj|), and (ii) the "plateau" near the maximum of 7 m (T) 
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FIG. 2: (color online), (a) Resistivity p m (T) (solid lines), and, (b) 
decoherence rate 7m(T) for 2S = n = 1, 2, 3; p m (0) = 2rp/nhi/Q, 
7°, = 2/-kFlvo, where p is the residual resistivity, vo the density of 
states per spin and channel, r the elastic scattering time, and 7™ ax 
is the maximum value of 7 m (T). We defined the Kondo scale Tjf 
for each S via p m {T^) = p m (0)/2. Dashed lines in (a) show that 
the empirical form p m (T)/p m (0) « /s(T/T|) with / s (a;) = (1 + 
(2 1/qs - l)x 2 )" as , used to fit experimental to NRG results for 
5 = 1/2 (H], also adequately fits the NRG results for S = 1 and 
S = 3/2. 



increases slightly in maximum height 7™ ax and significantly 
in width. These changes turn out to be sufficient to identify 
the proper value of S when comparing to experiments below. 

Comparison with experiment. — We compared our theoret- 
ical results for p m (T) and "fm(T) to measurements on quasi 
1 -dimensional, disordered wires, for two AgFe samples jl ill . 
(AgFe 2 and AgFe 3 having 27 ± 3 and 67.5 ± 7 ppm Fe im- 
purities in Ag, respectively), with a Kondo scale Tk ~ 5K 
(for S = 3/2, see below). These measurements extend up to 
T < Tk allowing the region T/Tk < 1 of the scaling curves 
in Fig. [2] to be compared to experiment. At T > Tk ~ 5K 
(i.e. T/Tk > 1) the large phonon contribution to the deco- 
herence rate prohibits reliable extraction of j m (T) for our Ag 
samples (see below) . In order to compare theory and experi- 
ment for temperatures T/Tk > 1, above the maximum in the 
decoherence rate, we therefore carried out new measurements 
on a sample (AuFe 3) with 7 ± 0.7 ppm Fe impurities in Au 
with a lower Kondo scale Tk ~ 1.3Kbut, as discussed above, 
described by the same Kondo model. Combining both sets of 
measurement thereby allows a large part of the scaling curves 
in Fig. |2]to be compared with experiment. 

Following ll ill , we subtract the electron-electron contribu- 
tion 124ll from the total resistivity p, yielding Ap due to mag- 
netic impurities (m) and phonons (ph): 
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FIG. 3: (color online). Measured resistivities Ap(T) (symbols) fit- 
ted to equation (|3j (lines), for n = 2S = 1,2 and 3, in the range 
below the onset of the phonon contribution, but above 100 — 200 mK 
l25ll . Specifically, we used 0.1 - 1.6 K for AuFe and 0.29 - 5.9 K 
for AgFe (arrows). The curves for AgFe 2 and AuFe 3 have been 
offset vertically by 0.25 and 0.75, respectively. The inset gives the 
Kondo scales for AgFe and AuFe extracted from the fits. Esti- 
mates of the unitary Kondo resistivities for n — 1,2 and 3 (in units 
of nO.cm/ppm) yield p m (0) = 0.041, 0.047 and 0.049 for AgFe 
(averaged over the two samples), and 0.23, 0.26 and 0.27 for AuFe, 
respectively. 



Here S is an (unknown) offset 026ll and Ap(T) is expressed per 
magnetic impurity. For temperatures low enough that p p h(T) 
can be neglected, Ap(T) — S corresponds to the theoreti- 
cal curve p m (T) - p m (0)f s (T/Ti) [cf. caption of Fig. 0, 
where p m (0) — Ap(0) — 5 is the unitary Kondo resistiv- 
ity. Fig. [3] illustrates how we extract the Kondo scale 
and Pm(0) from the experimental data, by fitting the Kondo- 
dominated part of Ap(T) in a fixed temperature range (spec- 
ified in the caption of Fig. |3]l to the NRG results of Fig. Ufa), 
using the Ansatz 



Ap(T) » 8 + (Ap(0) - S)f s (T/Tl) 



(3) 



Ap(T) = p m (T) + p ph (T) + S . 



(2) 



Such fits are made for each of the fully screened Kondo mod- 
els, using Tjf and 5 as fit parameters. Importantly, the val- 
ues for Tj| and p m (0) obtained from the fits, given in the 
inset and caption of Fig. [3] respectively, show a significant 
5-dependence: both Tjf and p m (0) increase with S, since the 
slope of the logarithmic Kondo increase of the theory curves 
for p m [cf. Fig. 12 decreases significantly in magnitude with S. 
Nevertheless, all three models fit the Kondo contribution very 
well, as shown in Fig. [3] so a determination of the appropriate 
model from resistivity data alone is not possible. 

To break this impasse, we exploit the remarkably sensis- 
tive iS-dependence of the spin-fiip-induced decoherence rate 
7 m (T). Fig. |4] shows the measured dimensionless decoher- 
ence rate 7m(T)/7™ ax for Ag and Au samples (symbols) as 
function of T/Tg for S = 1/2, 1 and 3/2, using the T| val- 
ues extracted from the resistivities, together with the corre- 
sponding parameter-free theoretical predictions (lines), taken 
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FIG. 4: (color online). Comparison of the measured (symbols) 
and theoretical (lines) results for the dimensionless decoherence rate 
7m(T)/7S ax as function of T/Tf, using 5 = 3/2. Insets show 
comparisons to S = 1/2 (a) and S = 1 (b). Tjf for AgFe and 
AuFe was extracted from the resistivities (inset of Fig. 0, while 
7™™ was determined as the average plateau height in the region 
T/T'^J 2 e [0.7, 1.35]. Typical error bars are shown for S = 3/2. 
They grow with increasing temperatures due to the increasing diffi- 
culty of subtracting the growing phonon contribution to the decoher- 
ence rate. 



from Fig. |2jb). The agreement between theory and experi- 
ment is poor for S = 1/2, better for S — 1 but excellent for 
S = 3/2, confirming the conclusion drawn above from ab 
initio calculations. The dependence on S is most strikingly 
revealed through the width of the plateau region (in units of 
T/Tg), which grows with S for the theory curves but shrinks 
with S for the experimental data (for which Tjf grows with 
S), with S — 3/2 giving the best agreement. 

Conclusions. — In this Letter we addressed one of the fun- 
damental unresolved questions of Kondo physics: that of de- 
riving and solving the effective low-energy Kondo model ap- 
propriate for a realistic description of Fe impurities in Au and 
Ag. Remarkably, for both Ag and Au samples, the use of a 
fully screened S — 3/2 three channel Kondo model allows a 
quantitatively consistent description of both the resistivity and 
decoherence rate with a single Tk (for each material). Our re- 
sults set a benchmark for the level of quantitative understand- 
ing attainable for the Kondo effect in real materials. 
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